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Abstract. The Ritt problem asks if there is an algorithm that tells whether 
one prime differential ideal is contained in another one if both are given by their 
characteristic sets. We give several equivalent formulations of this problem. In 
particular, we show that it is equivalent to testing if a differential polynomial 
is a zero divisor modulo a radical differential ideal. The technique used in the 
proof of equivalence yields algorithms for computing a canonical decomposition 
of a radical differential ideal into prime components and a canonical generating 
set of a radical differential ideal. Both proposed representations of a radical 
differential ideal are independent of the given set of generators and can be 
made independent of the ranking. 



1. Introduction 

The Ritt problem is an algebraic treatment of the following phenomenon occur- 
ring with differential equations. One can show that the solution set to the differ- 
ential equation y'^ — 4y = consists of a family of parabolae and the zero special 
solution. The solution set to the equation y'^ — Ay^ = consists of hyperbolae and 
the zero solution as well. The major difference between these two solution sets is 
that in the first case the zero solution is not a limit of the family while in the second 
case it is. In the algebraic language, the first equation does not generate a prime 
differential ideal, while the second one does. The Ritt problem is a generalisation 
of this phenomenon to the case of several non-linear algebraic PDEs. 

More precisely, consider a differential polynomial ring over a differential field of 
characteristic zero. J.F. Ritt has shown that every radical differential ideal can 
be uniquely represented as a finite intersection of its essential prime components, 
which are the minimal prime differential ideals containing it. We call the problem 
of computing the characteristic sets of the essential prime components, given a 
finite set of generators of a radical differential idea, the generalised Ritt problem. 
As of today, the problem remains largely unsolved, with some important special 
cases having been treated by J.F. Ritt, H. Levi, E.R. Kolchin, and R.M. Cohn 

[HIillIiHlElll]. 

We give a precise formulation of this problem as a computational problem, em- 
phasise its dependence on the differential field, and then reformulate the problem to 
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make it ficld-indcpendcnt. One of the versions of the Ritt problem that we diseuss 
here is given in terms of zero-divisors modulo an ideal. In particular, when we test 
primality by definition for a pair of polynomials each not in the ideal we need to 
determine whether their product is in the ideal. This sentence has two universal 
quantifiers, one for each polynomial. We show that it is sufficient to quantify with 
respect to just one polynomial. 

To illustrate these ideas we also exhibit an algorithm whose input is a finite set 
F of differential polynomials and the output is a collection of prime differential 
ideals given by their characteristic sets such that: 

• the radical differential ideal / generated by F is the intersection of the 
prime differential ideals and 

• this set of prime differential ideals is uniquely determined just by the ideal 
/ and does not depend on the chosen set of generators F. 

We hope that the approach we propose in this paper will bring us closer to deter- 
mining decidability of the Ritt problem. 

2. Basic definitions 

Let k be a differential field of characteristic zero with derivations A = {di, . . . , (5m}- 
Let Y — {yi, ■ ■ ■ ,yn} be a finite set of differential indeterminates and 

QY := {SY ■■■■■ Ufc ^ 0, 1 fc s; m; 1 s; J < n} . 

The ring of differential polynomials kjl"} is the ring of commutative polynomials 
k[0y] with the natural structure of a A- ring. 

Let / be a radical differential ideal in k{y}, that is, a radical ideal stable under 
the action of A. A prime differential ideal P is called an essential prime component 
of / if / C P and there is no prime differential ideal Q satisfying I Q Q C. P. 
According to [18j , every radical differential ideal / has finitely many essential prime 
components and equals their intersection; every prime differential ideal containing / 
must contain an essential prime component of /; and every set of prime differential 
ideals whose intersection equals / contains the set of essential prime components 
of / as a subset. The generalised Ritt problem is: given a radical differential ideal, 
find its essential prime components. 

To formulate it as a computational problem, we need to fix an encoding of the 
input (a radical differential ideal in k{F}), and of the output (a finite set of prime 
differential ideals in k{y}). A radical differential ideal is usually given by a finite 
generating set of differential polynomials. To be able to compute with them, we 
have to assume that k is a computable differential field. The following definition is 
directly derived from [ITl Definition 5], see [15] : 

Definition 1. A differential ring R is called computable, if there exists an embed- 
ding 

i: R^N 

such that 

(1) The image i{R) is a recursive subset ofN 

(2) The addition, multiplication, and derivation functions induced on this im- 
age by the addition, multiplication, and derivations on R are computable 
functions. 
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The embedding i in the above definition is caUed an admissible indexing of R. If 
k is a computable differential field, then the ring of differential polynomials k{y} 
is a computable differential ring: the value of an admissible indexing j : k{l"} N 
can be effectively computed for each differential polynomial / G k{y}, given the 
indices of its coefficients w.r.t. the indexing i : k ^ N 17i|- 

Having fixed such an indexing j, we can unambiguously say that a differential 
polynomial / € k{y} is given, meaning that actually j{f) is given. A prime 
differential ideal can be given either by a generating set (since it is radical), or by 
a characteristic set. The latter is defined as follows. Define a ranking on QY as a 
total order relation < satisfying: 

u < Su, u < V ^ Su < 6v 

for all u,v E QY, S E A. For a differential polynomial / e k{y} \ k, call the 
derivative of the highest rank effectively present in / the leader of /, denoted 
u = ld< /. If the degree of / w.r.t. u is d, then u'^ is called the rank of /, denoted 
rk< /. The ranking then extends to the set of ranks: 

< ui < U2 or (ui = U2 and di < c?2) 

and to the set of all finite sets of ranks: i?i < R2 iff either i?i — R2, or else the 
minimal element of the symmetric difference (i?i \ R2) U (i?2 \ Ri) belongs to Ri. 
All three relations: < on the set of derivatives, < on the set of ranks, and < on the 
set of finite sets of ranks, are well-orderings. 

Rankings are introduced, in order to be able to effectively represent prime differ- 
ential ideals by their characteristic sets, which allow to test membership and solve 
other algorithmic problems. The latter are defined as follows. 

A differential polynomial / is said to be reduced w.r.t. a differential polynomial 
g, if for any derivative operator 9 E 0, the degree of / w.r.t. to the the leader of 
9g is less than that of 6g. A set of differential polynomials is called autoreduced, if 
each element of this set is reduced w.r.t. the rest. Every autoreduced set is finite. 
Among all autoreduced subsets of a given set of differential polynomials X, choose 
one with the least set of ranks (which exists, because < is a well-ordering on the 
set of finite sets of ranks). This set is called a characteristic set of X. Note that 
the characteristic set may not be unique, but its set of ranks is. 

Consider a differential polynomial / as a polynomial in its leader u. Then the 
leading coefficient of / is called its initial, and the initial of any proper derivative 
9f, where 9 E 8\{1}, is called the separant of /. For a finite set C of differential 
polynomials, denote by He the product of the initials and separants of its elements. 
It is a well-known fact (see e.g. [12]) that every prime differential ideal P can be 
represented by its characteristic set C as 

that is, as the differential ideal generated by C and saturated by the initials and 
separants of C. We will extensively exploit this representation in the paper. 

3. Main result 

3.1. Different ways of posing the Ritt problem. We give a list of equivalent 
formulations of the Ritt problem below and will use one of them to attack the 
problem. We note that the equivalence takes place only for the fields k with a 
splitting algorithm, that is, an algorithm which, given a univariate polynomial / E 
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k[a:;], determines whether / is irreducible. This requirement is usually imposed 
in the context of polynomial factorisation problems over k (see, e.g., [20l [T6]). 
Examples of computable fields with splitting algorithm include the fields of rational 
numbers and rational functions, and their algebraic closures. An example of a 
computable field without splitting algorithm is given in [16j : take a computably 
enumerable but non-computable subset S of the natural numbers, and consider 
the field Q[y/p^ : n € S], where p„ is the n-th prime. By considering the field of 
fractions over this field, we obtain a non-trivial differential field without splitting 
algorithm. As in '15', Section 5], we suggest to think of the Ritt problem, especially 
of its third and fourth formulations in the theorem below, as a generalisation of the 
polynomial factorisation problem. 

Theorem 1. The following problems are equivalent over a computable differential 
field of characteristic zero with a splitting algorithm (an algorithmic solution for 
any one of them will also provide algorithms for the rest): 

(1) Given a characteristic set of a prime differential ideal, find a set of its 
generators. 

(2) Given the characteristic sets of two prime differential ideals Ii and I2, 
determine whether Ii C /2 - 

(3) Gompute a non-redundant prime decomposition of a radical differential 
ideal. 

(4) Given a radical differential ideal specified by a set of generators, determine 
whether it is prime. 

(5) Given a radical differential ideal specified by a set of generators, compute its 
prime decomposition together with the generators of the prime components 
as radical differential ideals. 

(6) Given a radical differential ideal I specified by a set of generators and a 
differential polynomial f , determine whether f is a zero-divisor modulo I . 

Proof. 1 2 Assume that we have an algorithm for finding generators of a prime 
differential ideal specified by its characteristic set. Applying this algorithm 
to the characteristic set of /i, compute its generator system Fi, i.e., Ii — 
{Fi}. Then h C h if and only if Fi C h. 
2 3 Given a radical differential ideal {F}, one can apply the Ritt-Kolchin algo- 
rithm [m Section IV.9], T^, Algorithm 5.5.15], pJi» Section 10] to compute 
its prime decomposition 

{F}^Pin...nPk, 

where each prime component is represented by its characteristic set. 

Note that the Ritt-Kolchin algorithm requires to test whether algebraic 
ideals of the form (C) : H^, where C is a coherent [121 Section III. 8] au- 
toreduced set, are prime and, if not, to find two polynomials not from the 
ideal, whose product belongs to the ideal. Apart from this test, the Ritt- 
Kolchin algorithm can be executed over any computable differential field of 
characteristic zero. 

The special case of C consisting of a single univariate polynomial shows 
that the existence of a splitting algorithm over k is necessary for the above 
primality test. We will show that it is also sufficient. For the proof of 
sufficiency of slightly stronger requirements see [131 Section 5.5]. 
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Indeed, if C is coherent autoreduced, then the ideal (C) : iJ^ is radical 
[21 [9]. A basis of this ideal can be computed over any computable field 
using Grobner bases [51 Section 4.4]. Then, the generators of the associated 
primes of the ideal can be computed over any computable field with the 
splitting algorithm [6]. If there is only one associated prime, the ideal is 
prime. Otherwise, by picking in each associated prime a polynomial that 
does not belong to the other associated primes (which can be done with 
Grobner bases over any computable field), we can find the required product. 

Assume that we have an algorithm for checking inclusion of prime differ- 
ential ideals. Then in the above decomposition we can remove all redundant 
components, i.e., those Pi that contain another Pj. The resulting compo- 
nents will be essential, since any prime decomposition of a radical ideal 
contains all essential components. Hence, we obtain the non-redundant 
(essential) decomposition of {F}. 

3 4 Assume that we have an algorithm for computing a non-redundant (=es- 

sential) prime decomposition of a radical differential ideal {F}. Then {F} 
is prime if and only if this decomposition consists of one component. 

4 1 Given a characteristic set C of a prime differential ideal P. Consider the 

algebraic ideals Ji = (C^*^) : H^, i — 0,1,2, . . ., where 

C« -{6*/ I ord6'<i, feC}. 

Let Fi be a system of generators of Ji (e.g., its Grobner basis). By the 
basis theorem, there exists an index i such that {Fi} = P. 

Assume that we have an algorithm for determining whether a radical dif- 
ferential ideal {F} is prime. Applying this algorithm to {Fi}, i = 0,1,2, .. ., 
find the least i such that {Fi} is prime. Then, since C is a characteristic 
set of {Fi}, we will have {Fi} — P, hence Fi generates P as a radical 
differential ideal. 

1^5 Given a radical differential ideal, one can apply the Ritt-Kolchin algorithm 
to compute its prime decomposition, where each prime component is rep- 
resented by its characteristic set. Assuming that we have an algorithm for 
computing generators of a prime differential ideal, given its characteristic 
set, we obtain the generators of the components. 

5 3 Assume that we have an algorithm for computing a prime decomposition 

of a radical differential ideal together with the generators of the prime 
components. Then we can check inclusion of the components as in the 
proof 1 2 and thus exclude redundant components as in 2 =J» 3. 
3 6 Assume that we have an algorithm for computing an essential prime de- 
composition / = Pi n . . . n Pfc of a radical differential ideal / — {F}. 

Lemma 1. A differential polynomial f is a zero- divisor modulo I if and 
only if it belongs to an essential prime component of I . 

Proof. Let / be a zero-divisor modulo /. Then by definition there exists a 
differential polynomial g ^ I such that fg e /. Since g ^ I, there exists 
an essential prime component Pi of / which does not contain g, yet fgGi 
implies /<? € Pi- Since Pi is prime, we obtain f £ Pi. 

Vice versa, let / be an element of an essential prime component of /. 
Let , . . . , Pi, be the essential prime components of / not containing /. 
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Then for the ideal I : f = {g \ gf E I}, we have 

/ : / = (Fi n . . . n no : / = (A : /) n . . . n (Pfe : /) = P,, n . . . n F,,. 

Since PiH . . .H Pk is a non-redundant prime decomposition of /, we obtain 
I : f ^ I, hence / is a zero-divisor modulo /. □ 

6^1 The proof of this implication is similar to [4 ^ 1]. 

Given a characteristic set C of a prime differential ideal P. Consider the 
algebraic ideals Ji = (C'*-*) : H^, i = 0, 1,2, . . .. Let Fi be a system of 
generators of J^. By the basis theorem, there exists an index i such that 
{Fi} = P. Note that the product He of initials and separants of C is not a 
zero-divisor modulo {Fi}, since 

P ^ {F^} C {F,} : C [C] ■.H^=P. 

Assume that we have an algorithm for determining whether a polynomial 
/ is a zero-divisor modulo a radical differential ideal {F}. Applying this 
algorithm to He and {Fi}, i = 0,1,2,..., find the smallest index i such 
that the He is not a zero-divisor modulo {Fi} (as we have shown above, 
such an index exists). Then {Fi} : H^ = {Fi}. On the other hand, 
{FJ : H^ ^P, since C C {FJ. Thus, {F,} = F. 

□ 

3.2. Canonical prime decomposition of a radical differential ideal. In this 
section we will show how one can produce a prime decomposition of a radical 
differential ideal / given by a finite set of generators F that does not depend on 
this particular choice of generators. A good candidate for this decomposition would 
be the essential prime decomposition but, again, it is still unknown whether the 
latter is computable. Our procedure is based on two simple observations contained 
in Lemma [2] and Proposition [TJ 

Assume that a ranking is fixed. Then a prime differential ideal has a canonical 
characteristic set uniquely determined by the ideal [TJ [HI [101 [H] ■ And given F C 
k{y} we can compute some prime decomposition 

k 

(1) I^{F} = f][C,]:H^, 

where Ci are the canonical characteristic sets of the corresponding prime differential 
ideals. 

Lemma 2. LetC he a characteristic set of the highest rank among Ci, . . . ,Ck- Then 
the ideal P := [C] : H^ is an essential prime component of I of the highest rank, 
among all essential prime components of F 

Proof. Let Q = [B] : Fg°, where B is the canonical characteristic set of Q, be an 
essential prime component of / such that P ^ Q. Because of this inclusion and the 
definition of the characteristic set as an autoreduced subset of the least rank, we 
have rkC < rkB. Moreover, since Q is an essential prime component of /, it must 
be among the ideals [Ci] : H^", i = 1, . . . , fc, which means that B = Ci for some /, 
1 ^ I ^ k. Thus, due to the choice of C, we have rkC > ikB. We conclude therefore 
that rkC = ikB. Hence, according to [51 Lemma 13], prime differential ideals F 
and Q are equal. □ 
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Proposition 1. Let I = {F} be a radical differential ideal, and let C = Ci, . . . ,Ci 
be a finite set of differential polynomials satisfying C C / C [C] : . Then 

/ = [C] : i/^ n / : Ci n . . . n / : Ci n {F U He}. 

Proof. Let 

f e[C]: ni -.CiH ...ni : Cin{FU He}. 

Then 

a-f£l, l^i^l 
and there exists h G H^^ such that 

h-fe[C]. 

We then have 

h-fef-[C]c{f-Cu...,f-Ci}cI. 
and, therefore, / G {F} : H^. Since due to, for example, [9l Proposition 2.1] 

{F} = {F}:H^n{FUHc}, 

we conclude that f E I. The other inclusion follows from / C [C] : H^. □ 

If in decomposition ^ there are several characteristic sets of the highest rank, 
say, Ci , . . . , Cq, and = , . . . , Ci^p, , 1 < i ^ then 

(2) / = [Ci] -.H^^nl: Cia n . . . n / : Ci,p, n {F U i7c J H . . . 

n [C,] : n / : n . . . n / : C,,p, n {F U i/cj . 

Note that Ci, . . . ,Cq and, therefore, all ideals in ^ are uniquely determined by /, 
that is, they do not depend on the choice of generators of /. Moreover, the ideals 

(3) / : Q,, 
and 

(4) {FUHeJ 

strictly contain /, 1 ^ i ^ q, 1 ^ j ^ pi. Indeed, all Cij are elements of essential 
prime components of /, therefore / : Cij 3 /. And, since H^ does not belong to 
the corresponding essential prime component [Ci] : H^ , it does not belong to /, 
and we have {F U He^} 2 I- 

By the Ritt-Raudenbush theorem, a strictly increasing chain of radical differen- 
tial ideals terminates. Therefore, by computing recursively the canonical, generator- 
independent prime decomposition of ^ and ^ , we obtain a generator-independent 
decomposition of the original ideal /. We note that the prime decomposition re- 
quired in step 11]) can be computed for the radical ideals given in a saturated form, 
as in ([3]). 



Admitting a slight abuse of notation in this formula, we denote by F U He the union of the 
set F and the singleton set containing He- 
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3.3. Ranking-independent canonical decomposition and generators. Even 
though the canonical prime decomposition computed by the method described in 
the previous section is independent of the given generators of the radical ideal /, it 
does depend on the choice of ranking. It is possible to obtain a ranking-independent 
canonical prime decomposition via the following modification. As in ([T]), compute 
any prime decomposition of /. Then, instead of extracting the characteristic sets 
of the highest rank w.r.t. a fixed chosen ranking (this step is ranking-dependent), 
take all prime components whose canonical characteristic sets have the highest rank 
w.r.t. some ranking, and consider all these characteristic sets instead of the above 
Ci, . . . ,Cq. This step can be accomplished by computing the universal characteristic 
set [7] for each prime component. The remaining steps are the same. 
Finally, from the canonical prime decomposition 

k 

1=1 

one can obtain a canonical set of generators of the radical differential ideal, which 
depends only on this ideal. For j = 0, 1, 2, . . . consider the algebraic ideal 

1=1 

and compute its Grobner basis Bj w.r.t. the lexicographic term order induced on 
power products of derivatives by the ranking (if a ranking-independent canonical 
set of generators is sought, compute the universal Grobner basis). Stop at the least 
value of j such that {Bj} — /, and output Bj. Such j exists, because / = Ujlo ^3- 
Note that the equality of two radical differential ideals given by generators can be 
checked by testing membership of generators of one ideal to the other ideal. 
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